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SUMMARY OF PITCH-DAMPING DERIVATIVES OF COMPLETE AIRPLANE 
ÀND MISSILE CONFIGURATIONS AS MEASURED IN FLIGHT 
AT TRANSONIC AND SUPERSONIC SPEEDS 


By Clarence L. Gillis and Rowe Chapman, Jr. 
SUMMARY 


Llongitudinal-damping data in the form of the pitching-moment damping 
derivatives (Сша + Cm) аге presented апа summarized from МАСА flight 


tests of rocket-propelled models and full-scale airplanes. The experi- 
mental data are compared with calculations and a discussion of each con- 
figuration is given. Detailed conclusions are precluded by the lack of 
systematic configuration changes and the relatively unknown effects of 
various factors. A general comparison of the results revealed the fol- 
lowing trends and general conclusions. Fhe contribution to the pitch- 
damping derivatives of wings having 45° or less sweepback is erratic in 

* the transonic region and may be either positive or negative, leading to 
possible dynamic instability of tailless configurations using such wings. 
No satisfactory method of calculating the damping of such: wings in the 

ы transonic region is available and а smoothly faired curve between sub- 
sonic and supersonic talculated values will probably be unconservative. 
Configurations having triangular and swept wings with approximately 
60° leading-edge sweep and rounded airfoil sections exhibited less vari- 
ation in damping in the transonic region. Calculated values of the pitch- 
damping derivatives for configurations with horizontal tails are more con- 
servative and give better agreement with experimental results if the 
distance from the trailing edge of the wing mean aerodynamic chord (instead 
of from the center of gravity) to the tail is used to calculate the down- 
wash leg effect. The damping caused by downwash effects arising from lift 
on a forward surface due to pitching velocity and rate of change of angle 
of attack may be appreciable if the forward surface is at a large distance 
from the center of gravity. 
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INTRODUCTION 


Successful accomplishment of aircraft flight at transonic and 
supersonic speeds has served to place increased emphasis on the dynamic 
behavior of aircraft. The type of behavior referred to here concerns 
motions of the aircraft as a whole as compared to motions of component 
parts usually referred to as flutter. Dynamic instability of the lat- 
eral and directional modes of motion is a problem that has occupied the 
designer's attention for some time. A trend toward deterioration of 
the damping of the longitudinal motion of aircraft operating in the 
transonic region has also been experienced. This decrease in damping, 
as measured by cycles to damp, occurs because of the high altitude at 
which such flights are usually conducted and because of an actual 
decrease in the dimensionless aerodynamic damping-in-pitch derivatives 
at transonic speeds for some aircraft configurations (refs. 1 and 2). 


The present paper summarizes the available experimental data 
obtained by the National Advisory Committee for Aeronautics on pitch- 
damping derivatives measured in flight at transonic and low supersonic 
speeds. Most of the data contained herein were obtained from rocket- 
propelled models of airplane or missile configurations flown at the 
Langley Pilotless Aircraft Research Station at Wallops Island, Va. 

Some data are also shown for several piloted &irplanes obtained from = 
flights at Ames Aeronautical Laboratory апа the NACA High-Speed Flight 
Research Station at Edwards, Calif. 


Some of the damping information presented herein has been published 
previously as parts of general longitudinel-stability investigations and 


some is, as yet, unpublished. In the present paper the damping data are _ 


compared with each other and with values calculated by the usual methods. 


SYMBOIS - 
Сү, lift coefficient, 1/98 : 
Cm pitching-moment coefficient, Pitching monent 
A aspect ratio = 
C constant in stability equation 
Ту mass moment of inertia about Y-axis, slug-ft^ 


Кү radius of gyration in pitch, Итън, ft 


тъ "9. т Я 
аа МлылА 22”. 
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L lift 

M Mach number 

S wing area (including area enclosed within fuselage), sq ft 

St horizontal-tail area (total included area), sq ft 

ү forward velocity, ft/sec 

W weight of body, lb 

b wing врап in equation for с, ft 

b damping constant in stability equations 

B local chord, ft 

b/2 

С mean &erodynamic chord, a | c? dy, ft 

g acceleration due to gravity, ft/sec? 

i tail length, measured from center of grevity of configura- 
tion to center of pressure of tail, ft 

m mass of body, W/g, slugs 

p free-stream static pressure, lb/sq ft 

q dynamic pressure when used in equations for Ст and Cm 
and equations (2) anā (3), i урм, 1Ъ/ва ft 

t time, sec 

Y reference axis through center of gravity of configuration 
perpendicular to plane of symmetry 

«й angle of attack, radians 

7 specific heat ratio (1.40) 


5 control-surface deflection, deg 
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angle of downwash, radians = е 


angle of pitch, radians 
rate of change of angle of attack, da/dt 
phase angie, radians я 2 


frequency of oscillation, radians/sec 


forward surface, based on its own area and chord 
rear surface, based on its own агеа and chord 
tail : - 


trimmed, or mean value 
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TEST AND ANALYSIS PROCEDURE 


Experimental data presented were obtained in free flight by the 
free-oscillation technique. Іп this test method, the aircraft is dis- 
turbed from a trimmed condition, usually by means of а rapid elevator 
deflection апа the resulting short-period oscillation is recorded ав 
the elevator is held fixed. The method of analysis of these oscilla- 
tions to obtain the static and dynamic stability derivatives is ade- 
quately covered in several references (such as ref. 3). For the 
present purpose, only that portion of the procedure dealing with the 
pitch-damping derivatives is of interest. 


Making the usual assumptions of constant velocity, level flight, 
and linear aerodynamic derivatives, a solution of the two-degrees-of- 
freedom equation of longitudinal motion of an aircraft can be obtained. 
For any appropriate quantity such as normal acceleration or angle of 
attack, the solution is of the form 


a = се? cos (ox + Ф) + ат (1) 


The constant b in equation (1) defines the damping of an oscillation 
and in terms of the aerodynamic derivatives is given by 


=- \2 
ви. Bor - (аа + ош} | A 


Equation (2) may be solved for the sum of the damping derivatives to 
give 


From the flight tests, therefore, the sum of the damping derivatives 
Cmq and Си, шау be determined if the damping constant b is measured 


and the lift-curve slope Cig is known. The damping factor b 13 gen- 


erally determined from the envelope of the curve defined by equation (1). 
This envelope is determined from the flight record of the appropriate 
measured quantity which, for rocket models, is angle of attack. Тһе 
equation for the damping factor is 


p = l98e Ao [Aa (h) 
to - tl 


р, 
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where дол and Лао аге the amplitudes measured from the mean value 

of о аб times tj and to. For rocket-model tests, normal accelera- 
tion and angle of attack were measured to provide the lift- curve 51орев. 
In some cases for the full-scale airplane tests reliable angle-of-attack 
information was not available from the flight test so the lifi-curve 
Slope was obtained from wind-tunnel tests and equations (1) and (1) | 
would be written in terms of the lift coefficient. Тһе test and analy- 
sis method described does not permit separation of the derivatives С 
and Сша. This is not а severe limitation, however, because the damping 


is always proportional to the sum of the two derivatives regardless of 
the flight conditions or mass characteristics. 


In free-flight tests, accurate measurement of the damping deriva- 
tive (Сша + Спа) is difficult for several reasons. First, as shown by 


equation (2), the total damping is composed of the damping-derivative 
term and the lift-curvé slope term, the relative magnitudes of which 
depend on the radius of gyration in pitch, The inaccuracies in 

Ста + Сша; 88 Obtaihed by solution of equation (2), are proportional to _ 


the relative contribution of the Ст term to the total demping, b. 
Аз ап example, the Cra. term contributed аз much ав two-thirds of the 


total damping in some of the rocket models, and in а full- Scale test of _ 
&n &irplane (ref. 1) the Ст term contributed about oné-half of the 


total damping. Present design trends indicate that the ‘proportion of 
damping contributed by the Ст term on future airplanes will more - 
nearly approach the proportion for the rocket models. Other factors 
which may contribute to inaccuracies in measuring Cm, + Cg, аге non- 
linear aerodynamic derivatives, disturbances due to gusts, and any other 
effects on the oscillation peaks which define the envelope of the curve. 
Gusts and nonlinear aerodynamic derivatives usually appear ав apparent | 
changes in the damping coefficients. 


CALCULATION METHOD 


The experimental damping derivatives presented herein are compared 
with calculated values from theoretical investigations wherever appli- 
cable, or experimental test data where such are available. It has fre- 
quently been assumed that all the damping on conventional airplane 
configurations is caused by the tail. А damping moment derivative С 
results from the additional angle of attack at the tail caused by the 
velocity of pitching about an axis through the center of gravity. A 
damping moment derivative Cm; results from the lag in downwash at the 
tail surface due to the finite time required for the downwash discharged 
at the wing to reach the tail surface (ref. 4). The damping derivatives 
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arising from these concepts аге given by the equation 


(даа + Спа), = - (2 + 3 (2 сы) (Ж) (5) 


where the downwash lag contribution is represented by the de/da factor 
in equation (5). А factor such as this assumes that the downwash lag 
term is characterized by downwash discharged at the center of gravity of 
the configuration. It appears that this downwash lag term might be more 
representative if it is assumed that the downwash is discharged at the 
trailing edge of the mean serodynamic chord of the wing rather than at 
the center of gravity. Modification of formula (5) in accordance with 
this concept gives 


(Cng * Ong), = -(2 + E SS) а (8) (6) 


where 1! is the length from the trailing edge of the mean aerodynamic 
chord to the center of pressure of the tail. 


For airplanes in which the tail contributes the largest part of the 
damping, equation (6) is satisfactory for an approximate calculation. 
It obviously fails for а tailless. airplane. For airplanes with swept 
Wings the damping contributed by the wing may be of appreciable magnitude 
at all speeds, and in the transonic region the wing damping may be of 
primary importance for wings of any plan form. Theoretical studies 
(refs. 5 to 8) and experimental data (refs. 9 and 10) show that, at 
transonic and low supersonic speeds, the wing itself may be dynamically 
unstable. Calculations of the damping derivatives should therefore 
include the wing even though the damping due to the tail may be the 
major factor. 


Additional increments in damping-moment coefficients arise because 
of the downwash on а rear lifting surface resulti from the lift on a 
forward surf&ce produced by the pitching velocity Ста) and the rate of 


change of angle of attack (94): The damping coefficients caused by 
this effect may be calculated ав follows: 


ACng = Tea д (т) 
гү 
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where 5 апа с are the quantities on which Cg is based. A similar 
expression сап be derived for ACms- Тһе sum of the two damping-moment 


coefficients is then: 


А (Cmq + Съ) = Un. Y = а За (Cra, + Ста) (9) 


As will be shown later this increment in the damping dérivatives is very 
small for conventional airplane configurations but сап be fairly large | 
for canard configurations. If the forward lifting surface is at а large 
distance from the center of gravity, the quantity a(g) can be obtained 
d 
гу 


to а close approximation by the expression 


The theoretical results used to calculate the pitch-damping deriva- 
tives for subsonic апа supersonic speeds are given in references 5 to 8, 
11, and 12. References 5 and 11 consider only untapered wings but the 
results were assumed to apply to tapered wings also. Wherever possible, 
wind-tunnel measurements of tail effectiveness апа downwash were used to 
calculate the damping caused by the tail. 


in calculating the derivative Cy; the simple downwash lag concept 
as given in equation (5) and modified in equation (6) was used at all 
speeds. Several refinements to these calculations have been considered 
in references 13 and 14 but for the present purpose these refinements 
were neglected because of the uncertainties in the experimental results 
caused by large possible experimental errors_in some cases and the rela- 
tively unknown but important effects of oscillation amplitude. 


RESULTS AND DISCUSSION 


The results of the flight measurements and the calculations are 
shown in figures 1 to 25 for airplane and missile configurations. А11 
information pertaining to one configuration is given in one figure; 
break lines on the drawings of the models indicate wedge and hexagonal 
airfoil sections. The details of the airfoil section for each model 
are given in table I. For consistency, damping coefficients in all fig- 
ures are based on the total area and mean aerodynamic chord of the wing | 
and are given for the rate terms (а and &) in radians per second. 

This method of presentation does not correspond to that с of the references 
from which the data were taken for some cases. 
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Scales to which the derivatives (Cmq + Cing, ) are plotted vary con- 


siderably between figures. It is important to note, as pointed out in 
reference 15, that a comparison of the absolute values of damping coef- 
ficients for different configurations has little significance. When 
presented, the curve designated “calculated (tail)" is the value obtained 
from equation (6), using wind-tunnel values of Clot and de/da. The 


calculated (tail) curve thus includes the downwash lag effect which is 

an indirect effect of the wing but acts on the tail. The direct damping 
of the wing, calculated from theoretical results for the conventional air- 
plane configurations was added to the calculated (tail) curve to give the 
curve designated "calculated (tail plus wing)." For the canard configu- 
rations, the calculated curve includes the damping increments given by 
equation (9) and the theoretical values of Cm, given by references 2 
апа б. There are insufficient data available to permit estimation of the 
effects of airfoil shape, oscillation amplitude, or oscillation frequency. 
Reference 16 shows some effects of oscillation amplitude and frequency on 
the damping derivatives (Сва + Ста.) Рог а 45° delta wing at subsonic speeds 
For reference, table II gives the reduced-frequency range for 211 models, 
and it may be noted that the frequency ranges for most of the data dis- 
cussed herein are below the ranges investigated in reference 16. 


Rocket-Propelled Models 


Model 1.- The experimental data for model 1 are contained in refer- 
ence 17 and, as shown in figure 1, exhibit a smooth variation through 
the transonic region. Theoretical values of Cmq + Coe and Cms, are 


also shown in figure 1 for the pertinent supersonic speeds. The positive 
value of Сту, at supersonic speeds decreases the Cm damping of the 
wing by approximately one-third. Although the experimental curve is con- 
siderably higher than the theoretical curve, the variations with Mach 
number are similar. The effect of the fuselage was investigated by the 
use of reference 18, which does not include the effect of the afterbody, 
and was found to be negligible. No data are available on the effects of 
the afterbody. 


Model 2.- Figure 2 shows the experimental дафа for two models of a 
tailless airplane configuration. The data for model 2 were taken from 
unpublished rocket-test results. No curve is faired since the data were 
at isolated Mach numbers and were obtained from sustained low-amplitude 
oscillations (Ла, = +0.259) at М = 0.91 and 1.24. The oscillation 
amplitude for the points at M = 1.34 was Лом t1.0". 


Model 2B (data previously unpublished) contained pulse rockets to 
disturb the model and, therefore, higher angles of attack were obtained 
in the test. The data obtained at higher angles of attack (Ла я ®1.0°) 
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confirmed the data for model А аб M = 1.34, which was for similar angle- 

of-attack variations. Model B also exhibited а continued oscillation of E 
low amplitude (50.259) similar to that which gave the positive damping | | 
coefficients for model A. Analysis of these low-amplitude oscillations, 
assuming time to damp to one-half amplitude 88 infinite, gave values | | = 
Of Сша + Cm, 88 varying РГгош 1.98 to 2.25. The values of damping 


obtained at higher amplitudes shown by the data points for model В and 

low emplitudes shown by positive data points for model A аге in reality 
boundaries for the damping coefficients where the coefficients for the `= 
system depends with certainty on the amplitude of the oscillation and 

possibly on the value of the mean angle of attack. С а 


In order to calculate the theoretical damping it was песезвагу to 
approximate the actual wing plan form by one more amenable to calcula- 
tion. Three approximations were tried: а swept tapered wing obtained "S 
by extending the leading and trailing edges to the root and tip, а: tri- 
angular wing having the same aspect ratio ав the actual wing, and а tri- 
angular wing having the same leading-edge sweep ав the actual wing. Ihe 
c&lculated curves are shown for 811 three approximations and best agree- 
ment was obtained between the calculated 52.50 triangular wing апа the 
flight data. 


Model 3.- Data are presented in figure 3 for two tailless models 
designed to have different wing flexibilities, These даңа аге from а 
program instituted to investigate the effects of wing flexibility on 
longitudinal stability. Тһе solid. line (model А) in figure 3 is experi- | 
mental data for а 9-percent-thick wing of Laminated wood and metal con- MEC 
struction having medium rigid characteristics and is 1аре1еа "least NE 
flexible wing." The dashed line is for а wing having the same geometric 
shape and thickness but with thinner metal inserts, hence, model B had а” 
more flexible wing and is labeled "flexible wing." The least flexible 
wing had а stiffness of approximately one-half that of а solid aluminum 
wing of the same thickness апа the flexible wing had а stiffness of | M 
approximately one-third that of a solid aluminum wing of the same thick- _ un 
ness. These stiffness ratios were determined by comparison of the influ- 
ence coefficients for the solid aluminum wing and the laminated wing. 


À calculation of the Ста damping of the wing was made.at М = 0.7 


by the method of reference 11. Тһе subsonic calculation gave the value Е 
of Сша ав approximately -16, which is far too high and unconservative — 


when compared with the experimental data. It is possible. that the disa- 
greement between calculation and experiment at subsonic speeds шау be u 
due to а positive Cms for the wing which is not accounted for in the 225 


computation. А computation of damping at М = 1.414 was made by using 
references 7 and 8 ара good agreement with the experimental data мав е 
obtained. 
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Model h.- Unpublished data for а 45° swept-wing tailless configu- 
ration are shown in figure 4. The wing was а 6-percent-thick section 
constructed of solid dural. Qualitatively, the shape of the curves for 
damping of models ЗА, 3B, and 4 show remarkable similarity for Mach num- 
bers from 0.90 to 1.15 with damping for both models exhibiting a charac- 
teristic drop in this region. The wing on model 1 was comparable in 
stiffness to the flexible wing of model 3. Data for model 4 were not 
obtained below a Mach number of 0.88; hence, it is not known if the damping 
decreased іп a manner similar to that for model 3. One calculated point 
for the wing alone is shown at a Mach number of 1.38. The calculated value 
is conservative but shows poor agreement with the experimental data. A 
calculation of body damping by the use of reference 19 showed the body con- 
tribution to be small at supersonic speeds because of the afterbody shape. 


Model 5.- The data for model 5 (fig. 5) are from reference 20. 
This model was a wingless fuselage-tail configuration that has been used 
for tests of а number of wing plan forms for supersonic airplanes 
(models 6 to 12). The damping coefficients in figure 5 are based on the 
wing area and chord of model 6 for comparison purposes. The theoretical 
values of Cmq + Сце agree fairly well with the experimental values, 


particularly at supersonic speeds, as does the calculated values of С 


of the tail surface (curve labeled "calculated Ста (tail)"). Тһе agree- 
ment of the theoretical and experimental values is somewhat fortuitous, 
however, because, as shown in figure 5, the theoretical values of C 

are considerably higher than those calculated from measured values of 
tail effectiveness. The theoretical values of Си; which are those for 
the isolated tail uninfluenced by downwash, are positive throughout the 


supersonic speed range considered and become rapidly larger as the Mach 
number decreases below about 1.2. 


Model 6.- The data for model 6 (fig. 6) are from reference 20. Two 
geometrically identical models were flown; model B having a steel wing 
(designated "rigid wing" in fig. 6) and model A, having an aluminum wing 
(designated "flexible wing" in fig. 6). АС transonic Mach numbers rather 
large changes in the experimental damping data occurred as contrasted to 
the smooth variation with Mach number of wing-off data (fig. 5) and calcu- 
lated tail contribution. Ав indicated by the calculated curve including 
the direct wing effects, part of the loss in damping near М = 1.0 can 
be explained by a large positive value of Ста + Сте, for the wing. This. 
large destabilizing effect of the wing is caused by the positive value 
of Cms, being much larger than the negative v&lue of . Another 


reason for the large damping changes near M = 1.0 сап probably be 
found in unaccounted-for downwash changes in this region. Іп calcu- 
lating the Сю portion of the tail contribution to damping it was 


assumed, in the absence of appropriate data, that the downwash variation 
With Mach number was similar to the variation of lift-curve slope, which 


1.2 
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did not show large changes (ref. 20). Since large changes in static 
stability for this configuration took place in the transonic region | 
(ref. 20), it is probable that the downwash variation with Mach number _ 
was considerably different from that assumed. | 


The increment in Сша + Crs, caused by the downw&sh resulting from 
the Ст and Ст, on the wing was calculated for model 6 and is shown 


in figure 6. This contribution to the damping is very small because the 
terms Ста and Cry, аге very small and, at least for supersonic speeds, 


are of opposite sign. Since this condition should be approximately the 

same for all configurations having the forward surface near the center 

of gravity, this contribution to the damping was neglected for all other 

such configurations herein, eS 


Model 7.- Model 7 was identical to model 6 except that the wing 

had an NACA 65( 06)А004.5 airfoil section instead of the sharp-edge = 
hexagonal section on model 6. Data for model 7 are shown in figure 7 
and were taken from unpublished data. Model Т, which had an aluminum 
Wing, showed less damping than does model 6A and a larger disagreement _ 
between the low-lift experimental curve and the calculated curve for the 
tail. The reason for the disagreement in the three sets of experimental 
data has not been established. 


Model 8.- Two experimental curves (fig. 8) were obtained for model 8 
(ref. 21): one at low lift coefficients and one at high lift coefficients. 
The calculated tail damping shows а large decrease at Mach numbers near 
0.93, as did the test data, and this is caused by the large decrease in = 
downwash at these Mach numbers (ref. 22). The direct effect of the wing 
was fairly small in the region where it could be calculated. The direct - 
contribution of the wing (as computed) is destabilizing at the supersonic 
speeds of the test. -It is perhaps significant that when the effect of 
the wing is included, the variations with Mach number of the experimental 
and calculated curves are, in general, similar at supersonic speeds. 


Model 9.- The damping data for two 45° swept-wing models are pre- . 
sented in figure 9 and are not previously published. Model A had a 
steel wing and is denoted as "rigid wing" in the figure. Model B had 
an aluminum wing and the curve is labeled "flexible wing" in figure 9. 
Calculated values of the damping shown for models A and B are based on 
experimental tail effectiveness. An additional curve labeled "calcu- 
lated (tail, theoretical)" is shown in the figure and is based on theo- 
retical tail effectiveness. The direct wing contribution to damping at 
pertinent supersonic Mach numbers was small and is not included in the 
calculated curves. It is noted that the difference in the experimental | 
damping curves for models A and B is in the opposite direction to that | 
difference in the calculated curves for the two models, which ів а E" 
center-of-gravity effect. А difference such ав this might possibly be 
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explained Бу a higher đe/da for the more rigid wing. Severe changes 

in downwash for the flexible wing might explain the high peak and rapid 
change that occurs between Mach numbers of 0.93 and 1.03 for the flexible 
wing. Peak values and severe changes such as these occur on a swept-wing 
airplane which is discussed in a subsequent section. 


Model 10.- The data for a 60° swept wing model are shown in fig- 
ure 10 as obtained from reference 23. The direct contribution of the 
wing to the total damping is illustrated by the increment between the 
total damping curve and the tail-damping curve (includes downwash lag). 
Two additional curves are shown in the figure to illustrate different 
techniques of computation. The curve labeled "calculated (tail, theo- 
retical)" uses theoretical tail effectiveness plus a downwash lag term 
from estimated downwash. The peak value of damping which occurs at 
М = 0.85 is not predicted by the calculations but is similar to a peak 
that occurs in the data of model 11. 


Model 11.- The experimental data presented in figure 11 for model 11 
are reported in reference 24. This modified triangular-wing model exhib- 
ited damping characteristics that showed а marked change between Mach 
numbers of 0.84 and 0.88 where the damping derivatives almost doubled in 
value in a manner similar to model 10. Calculated values for the tail 
and for the tail plus wing are also shown in the figure and show a smooth 
variation over the Mach number range. The difference in variation between 
the experimental and calculated curves and а comparison with the data for 
model 5 indicate that the rise in damping between M=0.80 and M = 0.90 
can be attributed directly to the wing or to increased downwash from the 


wing. 


Model 12.- The data presented in figure 12 for the 60° triangular- 
wing model are from reference 22. Two experimental curves are shown in 
the figure: one for the low-lift-coefficient range and the other for 
the high-lift-coefficient range. This difference in damping for high 
and low lift coefficients is primarily attributed to the nonlinear down- 
wash behind the triangular wing. The calculated (low lift, wing plus 
tail) curve shown is for dé/da evaluated at zero lift and for experi- 
mental tail effectiveness. Calculated damping coefficients, shown in 
figure 12, for the wing-plus-tail (high lift) configuration indicate 
that the increased damping at high lift can be primarily attributed to 
the increased downwash. 


Model 13.- Damping data for a supersonic airplane configuration are 
shown in figure 13. Experimental data are shown for two models; model A 
data were at high lift and model В data were at low lift. Curves for 
tail damping and total damping for model B are also shown in the figure. 
No computations were made for the high-lift model. The disagreement 
between the two sets of experimental data at subsonic speeds is believed 
to be due to the changes in downwash for high and low lift coefficients. 


уа»... 
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Good agreement was obtained between the computed and experimental data 
at low lift. From the experimental curve and the fact that it fairs 
into the calculated (tail plus wing) curve, it appears that the direct 
wing contribution to damping approaches zero at Mach numbers of approxi- 
metely 1.17 and 0.90 but generally does not give much destabilizing 
effect. It is noted that if the calculated (total damping) curve were 
extended to lower supersonic Mach numbers, the unstable contribution of 
wing Сп; would make the agreement between calculation and experiment 


worse. 


Model 14.- Data for a 45° swept-wing airplane model from reference 26 
are presented in figure 14. Calculated curves are shown for the tail 
damping using both the modified and unmodified tail length for the down- 
wash lag term. Flight data vary unpredictably through the transonic 
region and a comparison with the calculated curve for the tail indicates 
that the wing contributes both positive and negative damping over the 
region. А single data point (unpublished) for a similar tailless con- 
figuration shows the wing to have а positive (unstable) damping coeffi- - 
cient at М = 0.79. The calculated curve using the modified tail-length 
concept gives better average agreement with the experimental data and is 
more conservative than the unmodified formula. No total damping calcu- _ 
lated curve is shown since the wing contribution was small at the higher 
Mach numbers where it could be computed and could not be computed over 
the Mach number range where data were obtained. | 


Model 15.- The experimental data shown in figure 15 for model 15 
were taken from reference 2T. А curve of calculated tail damping shows 
fairly good agreement with the experimental data, both as to magnitude 
and general variation with Mach number. АС subsonic speeds the calcu- 
lated wing damping was appreciable, being about one-third of that con- 
tributed by the tail. For the supersonic speed range below а Mach number 
of 1.28 the calculated wing damping coefficient is positive апа the cal- 
culated damping of the complete configuration indicates & grester vari- 
ation with Mach number than is shown by the experimental data, 


Model 16.- Reference 28 and figure 16 contain the experimental data 
for model 16. Data are presented for models А апа B, identical models 
except for center-of-gravity location. The effect of the center-of- 
gravity location on the damping is small. Іп éontrast to most other | 
swept-wing configurations included herein, the variation with Mach num- 
ber of the experimental damping coefficients is smooth апа the damping 
increases ав the Mach number increases, The significance of these 
phenomena and whether they are related to the inverse taper аге not 
known. АС subsonic speeds, the calculated wing damping is about 50 per- 
cent of that calculated for the tail. Because of the inverse taper 
ratio on the wing, the methods of references Т and 8 are not applicable 


to this model. 
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Model 17.- The damping data for model 17, shown in figure 17, are 
from unpublished test results. This model had the same wing as the tail- 
less model 3 and had the same characteristic loss in damping near а Mach 
number of 1.0. Comparison of the experimental curve for model 17 with 
the calculated damping of the tail indicates that the wing contributed а 
positive (unstable) damping coefficient near M = 1.0 ав did model 3, 
although the quantitative values of wing damping for models 3 and 17 аге 
not directly comparable because of the different center-of-gravity loca- 
tions. Тһе agreement of the total calculated damping with the experi- 
mental values is very good аб subsonic speeds. The wing damping could 
not be calculated over the supersonic speed range of the tests. 


Model 18.- The results shown in figure 18 for а 609 trisnguler-wing 
missile configuration аге from reference 29. Damping data are presented 
for both the inline- and interdigitated-tail positions. As pointed out 
in reference 30 the marked difference in the characteristics for the 
inline and interdigitated tail сап be attributed to the difference in 
downwash at the tail by considering the rolled-up vortex sheet. The 
experimental curve for the inline-tail model has a variation with Mach 
number that is, in general, similar to that for the modified-triangular- 
wing model (model 11). A calculated curve for the tail is shown as well 
as a calculated curve for total damping coefficient (both curves being 
for the inline-tail configuration). Estimates of downwash and tail 
effectiveness based on experimental data were used and good average 
agreement with the flight data were obtained. A theoretical curve is 
shown and is generally nonconservative at the higher Mach numbers but 
shows good agreement with experimental data in the region near а Mach 
number of 1.1. | 


Model 19.- The damping data for а canard missile model shown in 
figure 19 are unpublished, The calcul&ted damping is considerably less 
than the experimental damping, but the variation with Mach number of the 
two curves is remarkably similar. The increments in damping coefficients 
caused by the downwash resulting from the canard lift effects СТ, 


and Сі, меге calculated and found to be appreciable as shown in fig- 
ure 19, This was true for all the canard configurations considered 
herein. The reason is that the value of Ста depends upon the distance 
of a lifting surface from the center of gravity (refs. 5 and 6) and, for 
the canard surface, becomes fairly large and of the same sign ав Cigs 
resulting in ап appreciable effect on the damping. 


Models 20 and 21.- The data for models 20 and 21 were taken from 
references 31 and 32 and are shown in figure 20. Models 20 and 21 were 
geometrically identical except that the cylindrical portion of the body 
between the wing ала the сапага was 10 inches longer for model 21 than 
for model 20. А11 the experimental curves in figure 20 show similar 
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variations with Mach number near а Mach number of 1.0, which differ con- 
siderably from the variations indicated by the calculated curves. For 
model 21 neither the theoretical damping nor that calculated from experi- 
mental data predicted the variation with Mach number of the measured | 
damping &t supersonic speeds. The calculated values for model 20 were _ 
considerably lower than the experimental values at supersonic speeds but 
showed a similar variation with Mach number. 


Model 22.- The experimental values for model 22 (fig. 21) were 
taken from unpublished results. Models А and B were identical except | 
for center-of-gravity location and the different lift-coefficient ranges 
covered by each model. Calculated values exhibit fair agreement with 
measured values only at supersonic speeds and low lift coefficients. 

The difference between the calculated values for models A and B was 
caused partly by the different center-of-gravity locations and partly 
by different values of downwash (shown by wind-tunnel data) over the two 
lift-coefficient ranges. Also shown in figure 21 are the values of 
damping contributed by the downwash resulting from the canard lift 
effects g^ ала Ста: ex баз 


Full-Scale Airplanes © B 


Airplane 1.- Damping data from flight tests of the .X-1 research air- 
plane are shown in figure 22. These dat& were obtained from reference 33, 
and reference 22 was used to obtain the calculated curves shown on the 
figure. The calculated damping agrees well with the experimental data 
both аз to magnitude and variation with Mach number. The calculated 
total damping curve is somewhat lower and agrees better with the experi- _ 
mental data than the calculated curve given in reference 33. This is 
probably caused by the use herein of the modified tail length for com- . 
puting the downwash lag effect (eq. (6)). Model 8, previously discussed, 
was equipped with the Х-1 wing. The low-lift-experimentel data for 
model 8 show a variation with Mach number ay eae to that for the 
flight-test data of the Х-1 airplane. и К - 


Airplane 2.- Тһе data for the X-4 airplane are shown in figure 23, 
ав Obtained from reference 2. Measured values are considerably lower 
than the calculated values and, for the highest altitude test, became 
zero at а Mach number of about 0.85. Аф а Mach number of 0.88, а sus- 
t&ined oscillation of small amplitude was obtained in the flight test 
which yielded a large positive value of Ста + Ста. Because this oscil- 


lation was of small amplitude, the value shown can be cofBidered valid 
only for small-amplitude motions. The measured values showed an effect 
of altitude, which, as stated in reference 2, may indicate that aero- 
elasticity has a large effect on the damping in this case, 


va № => 
vam dt zn 
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Airplane 3.- The results for the F-86 airplane in figure 24 were 
taken from reference 1. The sharp drop in the experimental curve at 
M= 0.92 was confirmed by a number of separate test points and corré- 
sponds to a similar irregular variation in the static stability 
parameter Спа. Definite reason for these fairly sharp changes was not 


established in the preliminary investigation of reference 1 but they are 
similar to those for other swept-wing configurations (models 3, 4, 9, 
апа 15). 


The estimated damping contributed by the tail and the tail plus 
Wing are shown in figure 24, Tail contributions were calculated from 
lift-curve slopes and downwash values obtained from transonic bump tests 
of 35° swept wings. Although reference 34 contains tail effectiveness 
data (0C,/dit) these were not used in calculating the damping because 
the values measured were about 25 percent lower than would be obtained 
by using the lift-curve slopes of an isolated 35° swept-tail surface. 
This is probably primarily caused by the fact that the inboard portion 
of the adjustable stabilizer (approximately 10 percent of the span) does 
not move. The lift contributed by the entire surface, however, is 
effective in producing damping. The damping contributed by the wing was 
appreciable at subsonic speeds for this configuration and a comparison 
of the calculated and experimental data indicates that the damping due 
to the wing is approximately zero at Mach numbers near 1.0. 


Airplane 4.- The damping data shown in figure 25 for the D-558-II air- 
plane were taken from reference 35. The calculated damping agrees fairly 
well with the experimental data at Mach numbers up to about 0.75, but as _ 
the Mach number increases, the experimental data show, first, a more rapid 
increase, and, then, a more rapid decrease than the calculated values. 

The more rapid decrease may be attributed, at least partially, to a proba- 
ble dynamically unstable contribution of the wing as shown by model 3 and 
airplane 2. к 


lt 


GENERALIZATION OF RESULTS 


A generalization of results on damping derivatives is limited by 
the number and types of configurations and the continuity of the experi- 
mental flight data available for summary. Some general trends are noted 
for each type of configuration апа а discussion is presented on the com- 
putation of damping. The general trends noted are limited for some of 
the cases by inconclusive evidence. 
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Tailless Configurations 
^ Four of the five tailless configurations (models 2,.3, and 4, and 
airplane 2) exhibited either very small negative values or positive 
values of Ста + Cms in the transonic region. The fifth configuration 


(model 1), with а 609 triangular wing, showed а fairly smooth variation 
through the transonic region with no indication of a decrease in damping 
Similar to that encountered on the swept-wing configurations. On model 2 
and airplane 2 the positive (unstable) damping coefficients were obtained 
from small-amplitude sustained oscillations and the evidence indicates 
that the damping is nonlinear with amplitude, being less аў the smaller 
amplitudes. The data shown for model 1 were obtained from fairly large- 
amplitude oscillations bút other models having the same wing on а slightly 
different fuselage have been flown with no pitch disturbance applied 
(unpublished) апа no sustained oscillations were encountered. 


Examination and comparison of the data for models ЗА, 3B, and 4 
show that the damping of swept wings is dependent on the flexibility. 
Also, these data show that similar swept wings exhibit similar varia- 
tions of damping with Mach number, characterized by a rapid decrease in 
damping above M = 0.85 which would not appear to be a critical Mach 
number effect. This same decrease in damping is evident: for airplane 2. 


Wingless Configurations B E 


Data are shown for only one wingless configuration (model 5) and 
indicate little variation in the damping coefficients with Mach number, 
It is reasonable to expect that the damping contributed by any lifting 
surface having а large moment arm from the center of gravity should vary 
with Mach number in а manner similar to its lift-curve slope if it is 
not operating in the flow field from a forward surface. о | 


Wing-Plus-Tail Configurations | — 


None of the configurations having tail surfaces exhibited any posi- 
tive values of Спа + Cms through the Mach number range tested, although 


several experienced large changes in the transonic region resulting in 
some fairly small negative damping coefficients. The contribution of an 
unswept wing to the damping derivatives of the total configuration is B 
usually erratic and пау be either positive or negative in the transonic 
region. Іп the majority of the cases the damping of the иизмеро пер: 
configurations was nonlinear with lift coefficient. 


In general the swept-wing configurations having 150 or less sweep 
also showed erratic changes in damping in the transonic region although 
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such changes were not as severe as those for the tailless swept-wing 
configurations because of the smaller contribution of the wing to the 
total damping. The one 60° swept-wing configuration discussed had a 
fairly smooth variation with Mach number near М = 1.0. The triangular- 
wing configurations generally had а somewhat smoother variation with Mach 
number than those with unswept or swept wings. Modifying or thickening 
the airfoil sections of the triangular wing (models 11 and 18) seems to 
result in somewhat less smooth variations in damping with Mach number. 


The results for model 17 show that the orientation in roll of the 
cruciform wings and tails with respect to each other had a significant 
effect on the damping. It has previously been shown in reference 36 
that the tail bending of a configuration similar to model 18 appreciably 
affected the damping in roll. It is probable that tail flexibility 
effects also influence the pitch damping. 


= 


Canard Configurations 


It appears that for the canard configurations shown, the experimental 
data are less amenable to accurate prediction, probably because of large 
interference effects and in the case of model 22 by nonlinearities in all 
the aerodynamic characteristics. Practically no experimental data are 
available at transonic speeds on the downwash and interference effects 
applicable to canard configurations. | 


The contribution to the damping coefficients of the downwash caused 
by the canard lift increments Cig апа Сү. was appreciable for all the 


canard configurations considered апа it appears that this effect should 
not be neglected for any configuration in which the forward lifting sur- 
Расе is at a large distance (either forward or rearward) of the center 

of gravity. | 


Agreement of Calculations and Experiment 


It has been generally shown that а computation of the damping con- 
tributed by the tail can be made and it will agree reasonably well with 
the experimental damping. If it is assumed that the downwash lag effect 
of the wing on the tail is that effect caused by the change in downwash 
in the finite interval of time required for air to travel from the 
trailing edge of the mean aerodynamic chord to the tail rather than from 
the center of gravity to the tail, then the damping due to the tail will 
exhibit better agreement with ЕЕ data and will be а more соп- 
servative estimate. 
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There exists по suitable means of calculating the wing contribution 
to damping in the Mach number region from 0.80 to 1.20 and a smoothly 
faired curve through this region from subsonic to supersonic speeds is 
of no value except possibly for triangular wings having rounded airfoil 
sections and which are not in the presence of а wake from а canard 
surface. 


The means of calculating the isolated surface contribution to 
damping by using references 5 to 8, 11, and 12 generally give values 
that are nonconservative but are in reasonable agreement with experiment 
if used within the limitations imposed by the theory. 


It should be noted that for computation of wing damping at super- 
sonic speeds, acceptable results are obtained only if : = Cmq and Спа 


are both used since the effects are generally in opposite directions. 
Computations of the damping derivatives for swept wings at lower super- 
sonic Mach numbers by use of references 7 and 8 are extremely limited by 
the theory. 


CONCLUSIONS 


From & study of the pitch-damping derivatives measured in flight on 
about 33 models and full-scale airplanes of widely differing configura- 
tion and geometric characteristics, the following general conclusions 
are offered. 

l. No suitable means exist for calculating the wing damping in the 
transonic region and a smoothly faired line between the subsonic and 
supersonic values will, in most cases, not conform to the experimental 
data and will probably give an unconservative result. 

2. The damping contribution of unswept and swept wings of 45° or 
less sweep is erratic in the transonic region and may be either positive | 
or negative, leading to possible dynamic instability for tailless 
configurations. | 


3. Configurations with triangular and swept wings having &pproxi- 
mately 609 sweep of the leading edge and rounded airfoil sections exhibit 
fairly smooth variations of damping in the transonic region. 


4. Calculations of the demping for configurations that have a hori- 
zontal tail surface appear to give the best agreement with the experi- 
mental data if the distance from the trailing edge of the mean aerodynamic 
chord of the wing to the tail is used for computing the downwash lag. 
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5. Damping moments resulting from the downwash on a rear lifting 
surface caused by the lift increments on a forward surface proportional 
to the pitching velocity and the rate of change of angle of attack may 
be appreciable if the forward surface is at а large distance from the 
center of gravity, ав on а canard configuration. 


It is believed that no detailed conclusions are warranted because 
(ОҒ the unsystematic nature of the data and the effects of relatively 
unknown factors such as experimental accuracy, nonlinearities with lift 


coefficient and oscillation amplitude, and the effects of oscillation 
frequency. 


Langley Aeronautical Laboratory, 


National Advisory Committee for Aeronautics, 
Langley Field, Va. 
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Figure 1.- Geometric characteristics and damping coefficients for a 
60° triangular wing aircraft configuration (model 1). 
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Wing 
А 2.01 
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Center-of-gravity location М.А.С., in. 21.90 
Model A | l6.7 percent M.A.C. 
Model B 9.9 percent М.А.С. 
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Figure 2.- Geometric characteristics and damping coefficients Рог а 
tailless aircraft configuration (model 2). 


30 


S, аа ft 

М.А.С., In Center-of—-gravíty location 
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Figure 3.- Geometric characteristics and damping coefficients for а 
150 sweptback wing апа body configuration (model 3). 
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Figure À.- Geometric characteristics and damping coefficients for а 
45° sweptback wing missile configuration (model 3), 
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Figure 5.- Geometric characteristics and damping coefficients for a 
wingless body and horizontal tail combination (model 5). 
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Figure 6.- Geometric characteristics and damping coefficients for an 
aircraft configuration having an unswept tapered wing with a 
. hexagonal airfoil section (model 6). 
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Figure 7.- Geometric characteristics and damping coefficients for 
an &ircraft configuration having an unswept wing with ап "E 
NACA 65( og A00h.5 airfoil section (model 7). 
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Figure 8.- Geometric characteristics апа damping coefficients for ап 
aircraft configuration having the X-1 wing (model 8). 


ЖЫР тағ» 


36 МАСА RM L52K20 


ДА 


Center-of—-gravity location 


Model А 24.8 percent М.А.С. = Wing 00 

Model 3 4.6 percent M.A.C. = ТА | 4.60 
S, sq ft 2.69 
М.А.С., іп. 10.05 


Horizontal Tail 
Same as model 5 


All dimensions are іп inches. жы 


Model А 
---------ұ----- Model В 


- С 
(са + ma) 


Caleulated (tail) VILI к 
“аде 
РЕ ||. _ 
BECP CO ES І” 
Сев ССС 


Caloulated (tail, theoretícal) à 


.7 .8 „9 1.0 1.1 1.2 1.5 1.4 1.5 В 


Figure 9.- Geometric characteristics апа damping coefficients for an . 
aircraft configuration having a 45° sweptback wing (model 9). 
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Figure 10.- Geometric char&cteristics and damping coefficients for an 
aircraft configuration having а 609 sweptback wing (model 10). 
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Figure 11.- Geometric characteristics and damping coefficients for ап. 
aircraft configuration having a modified triangular wing (model 11). 
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Figure 12.- Geometric characteristics апа damping coefficients for an 
&ircr&ft configuration having & 609 triangular wing (model 12). 


ы “. М t, > 
40 талы - ҠАСА RM 152К20 


bar 
І 


Wing 
Center-of-gravity location 20.5 А 5.0 
Model А O.4 percent M.A.C. 8.0 8, sq ft 4.25 
Model B -15.0 percent M.A.C. M.A.C., in. 15.2. 
"Span, in. 45.2 


Horizontal Tail 


A 5.0 
8, ва ft 0.79 
-М.А.С., іп. 6.54 


All dimensions are in inches. 


Model А 


TE aT NENET 
HNR 


10 


T „8 .9 1..0 1.3. 1.2 1.5 1.4 1.5 


Figure 13.- Geometric characteristics and damping coefficients for а 
supersonic aircraft configuration (model 13). 
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Figure 14.- Geometric characteristics and damping coefficients for an 
aircraft configuration having a 45° sweptback wing (model 14). 


42 NACA RM L22K20 


Wing 
А 4.00 
8, sq ft 5.56 
ма l = М.А.С., Іп. 14.64 
ЕЖ 
Center-of-gravity location 


10.9 percent М.А.С. ч Ще | =: 


7.00 Horizontal Tail 


A 5.72 
S, sq ft „958 
М.А.С., in, 6.27 


All dimensions sre in inches. 


e? «8 «9 2.0 1.1 1.2 2.5 1.4 1 5 
М 


Figure 15.- Geometric characteristics and damping coefficients of an 
aircraft configuration having а 40° sweptback wing (model 15). 
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Figure 16.- Geometric characteristics апа damping coefficients of an 
aircraft having а 409 sweptback wing that incorporates inverse 


taper (model 16). 


T МАСА RM L22K20 


Center-of-gravity location 


25,0 percent MeACe 


Horizontal tail 
A 4.00 
5, 34 ft (in chord 
p 


lane ) „90 
М.А.С., іп. 6.05 
All dimensions аге in inches. | 
Dihedral -209 


40 


МЕ Measured (low lift) 
реалов ag а) | | | 
Measured (Бі lift 
BINE (high э) 
SN 


ТІР Lr eene e 


ым ШОС 


Calculated (tail plus wing) 
d .8 .9 


1.9 2.1 2.2 2.5 


Figure 17.- Geometric characteristics and damping coefficients for an _ 


aircraft configuration having а 45° sweptback wing (model 17). 
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Figure 18.- Geometric characteristics апа damping coefficients of а 


609 triangular wing missile configuration (model 18). 
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Figure 19.- Geometric characteristics and damping coefficients for а 


missile configuration having & modified 609 triangular wing апа 
canard surface (model 19). | > 
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Figure 20.- Geometric characteristics and damping coefficients for two 
missile configurations having 60° triangular wings and canard surfaces 
(models 20 and 21). 


Wing 
Center-of-gravity location | K 3.02 
Model A Station 58.6 5 8. за ft 1.422 
Model B Station 40.4 Й.А.б., ine 9.20 
ids Canard 
A. 3.24 
8, ва ft „592 
М.А.С., іп. 5.55 


All dimensions are in inches. 


200 
ЕЕ Т sese à to aseo 
SE е. Calculated 

ГА ШЕ Experimental ЕШ 
T EA 4 NIMM 
в == аў AS 

Calculated 
TSR, Yap ER 


Increment due to canara Сі and A (1ow YS 
(high lif EAS 


"i .З .9 1.0 1.1 1.2 1.5 1.4 1.5 1.6 


40 


Figure 21.- Geometric characteristics апа damping coefficients for 
а missile configuration having essentially unswept wings &nd 
canard surfaces (model 22). 
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Х-1 airplane (airplane 1). 
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Figure 23.- Geometric characteristics апа damping coefficients for the 


X-4 airplane (airplane 2). 
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Figure 24,- Geometric characteristics and damping coefficients for the 
F-86 airplane (airplane 3). 
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D-558-II airplane (airplane 4), 
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